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Abstract 

The liquid-gas phase transition in strange hadronic matter is reexamined by using the new 
parameters about the A — A interaction deduced from recent observation of double hyper- 

nucleus. The extended Furnstahl-Serot-Tang model with nucleons and hyperons is utilized. The 
binodal surface, the limit pressure, the entropy, the specific heat capacity and the Caloric curves are 
addressed. We find that the liquid-gas phase transition can occur more easily in strange hadronic 
matter with weak Y-Y interaction than that of the strong Y-Y interaction. 
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I. INTRODUCTION 



Strangeness opens a new dimension to nuclear physics. Nuclear system with strangeness 
has many implications in astrophysics and cosmology and thus arouses much research inter- 
est. Nuclear systems with strangeness can be categorized into two species: strange quark 
matter (or strangelet) and strange hadronic matter(SHM). Quantum chromodynamics pre- 
dicts that the strangenss-rich system could be strange quark matter or strangelet, consisting 
of u, d and s quarks. Strangelet has been a hot research topic for both theoretical physics 
and experimental physics, but there is still no convincing experimental result to confirm 
its existence. The other possible existence of strangeness is SHM, in which quarks are lo- 
calized within nucleons and hyperons. The research of SHM is attractive because of the 
complexity of its structure and ingredients, especially, the phase transition including the 
quark deconfinement phase transition and liquid-gas(L-G) phase transition. 

This paper evolves from an attempt to study the L-G phase transition in SHM at finite 
temperature. In a previous paper we have studied the possibility of an L-G phase transi- 
tion in SHM utilizing the extended Furnstahl-Serot-Tang(FST) model. After investigating 
the mechanical and chemical stabilities of SHM, we found a critical pressure on the binodal 
surface for T = 13MeV and a limit pressure on the binodal surface for lower temperature. 
After examining the continuities of entropy and specific heat capacity, we found the L-G 



phase transition in SHM is of second order. Employin; 
model, the same conclusion has also been drawn in ref.|2 



;he chiral SU{3) quark mean field 



Recently, Takahashi et a/.P observed the A — A energy in a \/^He double hypernu- 



cleus. Their observation deduced the A — A energy to be AB 



AA = 1.01 ± 0.20t[j:}?Me1/, 
~ 5MeV from earlier 



which is much smaller than the previous estimation Ai?AA 
experiments 

mm 

. This leads to a significant decrease of potential Va^'* from about 
20MeV to 5MeV The new data imply that the hyperon-hyperon interact ion (Y-Y 

interaction) can be weaker than known before We employed two of many candidates 
of phenomenological SHM models 



model 



IG 
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Q, H Q, y , Q H Q 



, namely, extended FST 

Q 

12| and modified quark meson MC) model|9(, and the weak Y-Y 

interaction to reexamine the stability of SHM. We came to a conclusion that, while the 
system with the strong Y-Y interaction and in a quite large strangeness fraction region is 



more deeply bound than the ordinary nuclear matter due to the opening of new degrees of 



2 



freedom, the system with weak Y-Y interaction is rather loosely bound compared to the 

n 

latter [9]. This conclusion is not mo del- dependent on the extended FST model or MQMC 
model and true for both zero and finite temperature Q]. 

Noting that the stability affects on the L-G phase transition directly, the stability of SHM 
with weak Y-Y interaction is quite different from that with strong Y-Y interaction. It is of 
interest to investigate the influence of this difference on the L-G phase transition in SHM. 
To employ the extended FST model and reexamine the L-G phase transition in SHM with 
weak Y-Y interaction is the main purpose of this paper. We will show that there is still 
a second order L-G phase transition in SHM with weak Y-Y interaction, but the binodal 
surface, limit pressure and the equation of state are very different from those of strong Y-Y 
interaction. 

The paper is organized as follows. In Sec. II, the extended FST model is briefly introduced 
and the parameters for strong and weak Y-Y interactions are given. In Sec. HI we reexamine 
the L-G phase transition in SHM with weak Y-Y interaction. The last section is the summary 
of the main results. 



II. THE EXTENDED FST MODEL 



In refs. 
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12l |. we extended the original FST model to include not orily nucleons and 
0", uj mesons, but also A, S hyperons and strange mesons a* and (p. In ref.jl|, we used this 
model to stud y th e L-G phase transition in SHM. The details of this model can be found in 
Hereafter, we only give brief formulae, which are necessary for discussing 



lei to stud y th 



refs.j 

the L-G phase transition. 

In mean-field approximation, the Lagrangian of the extended FST model is presented as 
follows: 



Cmft = '4'N{h^d^ - g^^Nl Vq - Mn + gsNcro)'ipN 
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Qij are the coupling constants of baryon j to meson i field. 

By using the standard technique of statistical mechanics, the thermodynamic potential 
Q is obtained 

-2A;BT{^ln [1 + e-^^^*^*^)-'^^] + E In [1 + e-^^^tik)+^')\} (3) 
where (3 is the inverse temperature and V is the volume of the system 

El{k) = ^Mf + P (3) 
The effective masses of the hyperons and nucleons are 

M*^Mi- g,iao - ganc^l {i = A, 5) , (4) 
M*^Mi- gsiao (i = TV) (5) 

The mean-field values 00, Vq , ctq and ctq are determined by the corresponding extreme 
conditions of the thermodynamic potential. 
The baryon densities pBi is given by 

PBi = (il^t^i) = ^ j d^^"^ K (k) - fii {k)] (6) 

where — 2 ior i — N or E, — 1 ior i — A. The baryon and anti-baryon distributions 
are, respectively, expressed as 

n,{k) = {exp[(3{E:{k)-u,)] + l}-' (7) 

and 

riiik) = {expl/3{E*{k) + u^)] + 1}-^ (8) 
i/i are related to chemical potential /li by 

i^A = + gojAVo + g<t>A<Po, 

Us = i^s + S-o^hV^ + ff^S^O- (9) 



Since the system has equal number of protons and neutrons, and equal number of and 
S^, the chemical equilibrium conditions for the reactions A + A ^ ra + S" and A + A ^ 
are unified as 

2/Ua = /iiV + yUH- (10) 

Eq.(lO) implies that only two components of A^, A and S are independent. The SHM 
is regarded as a two-component system, for instance, the nucleon component and the S 
component. The strangeness fraction/5 is introduced as 

Pbk + '^Pbe 



Pb 
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(11) 
in the 



21 



22 



which plays the similar role as that of the asymmetric parameter a = {p. 
asymmetric nuclear matter. We can use the same method as that of refs. 
to address the L-G phase transition. 

Following the usual procedure of statistical physics, we can calculate the other ther- 
modynamic quantities from thermodynamic potential Q. For example, the pressure and 
entropy density are calculated by formulas p = —Q/V and S/V = —{dQ/d{l/P))v,^jV = 
{dp/d{l/(3)) 

v,tj.i, (70, Vo,<f>o,'^o- T^^^ results are expressed as below 



P 



V — 

4" 67r2 



dk 



Em 



[ni{k) + ni{k)]- Hg 



1-^ 



5*0 



0^0 



So 



'0 ' 



(12) 



k^ \ P\E: (k) - u,) exp [/? (Etik) - u,)] 



E*ik) \ [exp[/3iE*{k)-u,)] + l 



|2 



d\E:{k) + u,)exp[f3{E:{k) + u,)] 1 

[exp[/?(i?;(fc) + z.,)] + l]2 j K ) 

In the previous calculation Q], we employed the parameter set Tl given by ref. ^l|, and 
the earlier data of strong Y-Y interaction. In refs. 0, Q| , we have the new data of weak Y-Y 
interaction and adjusted the parameters from g'^*f^ = 48.31 and (7^,= = 154.62 for strong 
Y-Y interaction, to (7^,^ = 28.73 and (7^*^ = 129.06 for weak Y-Y interaction. The details 
of the parameters are summarized in Table 1. We will use the parameters for weak Y-Y 
interaction to address the L-G phase transition in the next section. 
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2.70 99.3 154.5 37.32 9.99 


9^* A (Strong 


;) (Strong 


r) 5^.A(Weak) r7^* = (Weak) 




48.31 


154.62 


28.73 129.06 





TABLE I: The parameters used in calculation 
III. L-G PHASE TRANSITION IN SHM WITH WEAK Y-Y INTERACTION 

In this section, we employ the extended FST model to investigate the L-G phase transition 
in the SHM with weak Y-Y interaction. 

Using the formulae in Sec. II and the parameters of Table 1, we calculate the pressure and 
study the equation of state first. The results for pressure vs. baryon density pb curves at 
T = lOMeV are shown in Fig.l, with weak and strong Y-Y interactions, respectively. We 
can see that the equation of state for two cases are quite different in two aspects. Firstly, 
with weak Y-Y interaction, when fs > 0.7, the pressure increases with density monotonically 
and there are no mechanical unstable regions. While with strong Y-Y interaction, there are 
always mechanical unstable regions for each fs- Secondly, with weak Y-Y interaction, the 
minimum of the pressure-density curve rises monotonically with fs, till fs > 0.7, where 
there is no minimum. While with strong Y-Y interaction, the minimum drops with fs and 
reaches the lowest value on the curve of fs = 1.3. When fs > 1.3, the pressure-density 
curves ascend with the increase of fs- These features are consistent with the behavior of 



3 and Q. 



energy-density curves and the stability of SHM in ref. 

To discuss the chemical instability of SHM with weak Y-Y interaction, we show the 
chemical potential isobars for nucleons and H against fs at temperature T = lOMeV" in 
Fig.2 for p = 0.05, 0.10, 0.20, 0.40, 0.50 and OmMeV fm-^ respectively. An inflection point 
could be spotted on the curve with the pressure = 0.50MeV fm~^. There are chemical 
unstable regions on curves with p < p*. This behavior is similar to that of strong Y-Y 
interaction but the inflection points are different for these two cases. 

The chemical, thermal and mechanical Gibbs equilibrium conditions for two separate 
phases require 

= T° = T (14) 
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/ij (t, p\ f^) = /ij^ (t, pG, /f ) , (g = AT, S), 



(15) 
(16) 



where the superscripts L, G denote the hquid and gas phases, respectively. We follow the 



procedures used in refs.j 



2l| to solve out Eqs.(14),(15),(16). Collecting up all 



the solutions and we get the phase boundary, i.e., binodal surface. The binodal surface at 
T = lOMeV are shown in Fig.3, where the dotted line refers to weak Y-Y interaction and 
the solid line to strong Y-Y interaction. We see from Fig.3 that the differences between 
two binodal surfaces are very remarkable. Though the limit pressures, above which the L-G 
phase transition cannot take place, exist for both two cases, the values are very different. 
The limit pressure for weak Y-Y interaction = 0.50MeV fm~^ is much higher than that 
of the strong Y-Y interaction as pf^^ = O.OdbMeVfm'^. The enclosed region of binodal 
surface with weak Y-Y interaction is much larger than that of the strong Y-Y interaction. A 
higher limit pressure and a larger region of phase separation imply that the system is likely 
to be less stable and the L-G phase transition can take place more easily. This result is of 
course very reasonable if we notice the conclusion given by refsjof Q| that the SHM with 
weak Y-Y interaction is less stable. 

To determine the order of the L-G phase transition, we examine the entropy density and 
the specific heat capacity of SHM with weak Y-Y interaction. By Ehrenfest's definition, the 
first order phase transition is characterized by the discontinuities of the first order derivatives 
of the chemical potential, such as the discontinuities of entropy and volume, while the second 
order phase transition unfolds the discontinuous behavior for the second order derivatives 
of the chemical potential, such as the specific heat capacity. The entropy per baryon is 



where S(T,p, fs) can be calculated from Eq.(13). 
The specific heat capacity is 



PB 



(17) 



as 



_9i_ 
67r2 



dk 



Em 



d'^Ui d'^Uj 
+ 



dT^ dT"^ 



where 
d rij 



{E*{k) - z/,)2(e^(^>*(^)~"0 - 1) - 2T{E*{k) - z/,)(e^(^«*('=)-"0 + 1) 



dT^ 



(18) 



(19) 



,2_ (,'^{Enk)+>y.) 



d Hi 
'dT^ 



{E*{k) + z/i)2(e^(^'(^)+"0 - 1) - 2T{E*{k) + z/i)(e'^(^'*(^)+"0 + 1) 



(20) 



The entropy vs. temperature curves for fs = 0.0,0.1 and 0.4 at p = O.llMeV fm'^ are 
shown in Fig. 4. We see from Fig.4 that the entropy evolves continuously through the phase 
transition when fs > 0.0. But when fs = 0.0, our model reduces to the one-component 
symmetric nuclear matter. The entropy becomes discontinuous. The L-G phase transition 
becomes the first order. This result is consistent with theoretical and experimental results 
about L-G phase transition in symmetric nuclear matter. 

The specific heat capacity vs. temperature curves for fs = 0.1 and 0.4 at p = 
0.11Mel^/m~^ are shown in Fig.5. The specific heat capacity curves are discontinuous 
for fs = 0.1 and 0.4. The continuity of entropy together with the discontinuity of specific 
heat capacity demonstrates that the L-G phase transition in SHM with weak Y-Y interaction 
is still of second order, as expected. 

One of the important signals for L-G phase transition is the investigation of Caloric 



curves 



23 



24 



25 



26l |. The Caloric curve reflects the relationship between excited energy of 



hadronic matter and the temperature. The energy of SHM is 
^ = E4 / dkeE*{k)[n,{k) + n,ik)]+Hg 
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So 



+gu;NVoPBN + {guAVo + 5'</.A0o)PBA + (S't^S^o + g^S(po)PBS 



(21) 



and the excited energy E* / pB is[ 

Xe\p\fs\T) + {l-\)e^{p^Js^,T) 



EVP 



B 



Ap^ + (1 - \)p^ 



(22) 



/ T=0 

We present the Caloric curves at pressure p = O.llMeV fm~'^ for fs = 0.0,0.1 and 0.4 
in Fig. 6. There are slopes during the L-G phase transition for fs = 0.1 and 0.4, while for 
fs = 0.0 there is a platform. This platform-shaped Caloric curve has been got in experiments 



25|. For fs > 0.0, the 



of the L-G phase transition in normal symmetric nuclear matter [ 
slopes mean that temperature does not remain constant during phase transition. This is 
the v ery feature of second order phase transition as was first pointed out by Miiller and 
Serot 
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To exhibit the behavior of Caloric curves of SHM with strong and weak Y-Y interactions 
clearly, we plot the temperature T vs. E* / pB curves at p = 0.11MeVfm~^, fs — 0.4 
for strong and weak Y-Y interactions in Fig. 7. We find the temperature region of phase 
transition for weak Y-Y interaction is lower than that of strong Y-Y interaction and the 
slope for weak Y-Y interaction is more acclivitous. Perhaps these distinctions can provide 
us a possibility to check the Y-Y interaction be strong or weak, and then the A — A energy. 

IV. SUMMARY 

We have used the new parameters about the weak Y-Y interaction derived from recent 
experiment data to reexamine the L-G phase transition in SHM. We have employed the 
extended FST model with nucleons and A, S hyperons and simplified the SHM into a two- 
component system to discuss its L-G phase transition. We come to the following conclusions: 

1. No matter the Y-Y interaction is strong or weak, the L-G phase transition can take 
place in SHM. The L-G phase transition is of second order for /s 7^ 0, the entropy 
is continuous and the heate capacity is discontinuous at the phase transition point. 
But for fs — the L-G phase transition becomes the first order. The entropy is 
discontinuous. 

2. Due to the difference between the stability of SHM with strong and weak Y-Y inter- 
actions, the physical features of L-G phase transition, including the limit pressure, 
the regions enclosed by the binodal surface, the Caloric curves, the slope of phase 
transition curves are very different for weak and strong Y-Y interactions. 

3. After reexaming the features of SHM, we come to a conclusion that the L-G phase 
transition can occur in SHM with weak Y-Y interaction more easily than that with 
strong Y-Y interaction. 
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FIG. 1: Pressure as a function of baryon density at temperature T = lOMeV for various strangeness 
fractions fs with weak and strong Y-Y interactions. Note that the scales of the two figures are not 
the same. 
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FIG. 2: Chemical potential isobars as functions of fs at T = lOMeV. The curves have pressures 
p = 0.05,0.10,0.20,0.40,0.50 and OmMeVfm-^ respectively. The inflection point(A) is at = 
0.50MeVfm-^. 
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FIG. 3: Binodal surface at T = lOMeV. The dotted hne is the biondal surface for weak Y-Y 
interaction, while the solid line is that for strong Y-Y interaction. The binodal surfaces are cut off 
at limit pressure pum = 0.50MeVfm~^ and pum = 0.095MeV fm~^ respectively. 
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FIG. 4: Entropy as a function of temperature at constant pressure p = 0.11MeVfm~^ for 
strangeness fraction fs = 0.0,0.1 and 0.4. For fs = 0.1 and 0.4, tlic entropy evolves continu- 
ously through the phase transition, while for fs = 0.0 the entropy is discontinuous. 
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FIG. 5: Specific heat capacity as a function of temperature at constant pressure p = 0.11MeVfm~^ 
for strangeness fraction = 0.1 and 0.4. For fs = 0.1 and 0.4, the entropy is discontinuous through 
the phase transition. 
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FIG. 6: Temperature as a function of excited energy at constant pressure p = 0.11MeVfm~^ for 
strangeness fraction fs = 0.0,0.1 and 0.4. For fs = 0.1 and 0.4, there are slopes during the phase 
transition, while for fs = 0.0 there is a platform. 
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FIG. 7: Temperature as a function of excited energy at constant pressure p = O.llMeVfm ^ for 
strangeness fraction fs = 0.4. 
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